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WASHBURN OBSERVATORY 

UNIVERSITY OF WISCONSIN-MADISON 

475 NORTH CHARTER STREET 
MADISON, WISCONSIN 53706 
TELEPHONE: (608)262-3071 

June 11, 1987 

Dr. Yoji Kondo 
Technical Officer 

NASA-Goddard Space Flight Center 
Greenbelt, Maryland 20771 

RE: NAG 5-832 Semi-Annual Progress Report 


Dear Dr. Kondo, 


This letter is intended as a report of the progress made on the NAS A-IUE project title A Study of 
the Fundamental Characteristics of 2175A Extinction with J. A. Cardelli and B. D. Savage 
as principle and co-investigators respectively. The data for this project were obtained with the 
International Ultraviolet Explorer satellite from the control center at NASA-Goddard Space Flight 
Center in Greenbelt, Maryland over two days in July and again in September of 1986. The principle 
investigator travelled to Goddard on both occassions. 


The purpose of the project is to study the characteristics of interstellar extinction in the region of the 
2175A feature for lines of sight which appear to exhibit unusually weak ultraviolet extinction. The 
objects that have been observed appear in Table 1. The analysis is based upon a parameterization of 
the observed extinction via fitting specfic mathematical functions in order to determine the position and 
width of the 2175A feature. The data are currently being analyzed and a publishable manuscript is in 
progress. The results so far indicate that even for lines of sight which exhibit unusually weak 2175A 
extinction, the central position is remarkably constant being generally within the limits that have been 
found for a large number of other lines of sight showing a wide range of observed extinction 
(Fitzpatrick and Massa 1986, Ap. J., 307, 286.). The width of the 2175A features observed for our 
objects, although generally within the range of variance of the sample in the above reference, tend to 
be slightly larger on average than the standard feature width. This seems to indicate that features 
associated with dark interstellar clouds may be broader on average than those associate with low 
density interstellar matter. 
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For one of our objects, HD 62542, the 2175A feature appears to be one of the weakest yet 
observed for any line of sight. The line of sight intersects a cloud of interstellar matter that appears to 
be exposed to both dense conditions as well as strong stellar winds. The analysis of the data for this 
line of sight have proven to be extremely interesting and a paper on this subject titled Two Lines of 
Sight with Exceedingly Anomalous UV Interstellar Extinction has already been submitted 
to the Astrophysical Journal. A copy of the manuscript has been enclosed. The analysis of the data 
has shown that the central position of the UV extinction feature is shifted from the standard value of 
2175A to 21 10A which represents the largest observed shift yet observed. The unusual physiology of 
this line of sight appears to have produced one of the most dramatic examples of the modification of 
the 2 175 A feature carrier yet observed. 


TABLE 1 PROGRAM OBJECTS 


Name 

V 

Spectral Tvpe 

EfB-Vl 

ABump 

HD 45314 

6.64 

09e 

0.46 

-1.98 

HD 149452 

9.05 

09 V 

0.90 

-1.62 

HD 191396 

8.13 

B0.5 II 

0.53 

-2.34 

HD 248893 

9.69 

B0 II/III 

0.74 

-2.02 

HD 199021 

8.45 

BOV 

0.86 

-1.87 

HD 33203 

6.02 

B2II 

0.96 

-3.18 

HD 164073 

8.04 

B3 ffl/TV 

0.24 

-3.16 

BD+45 3360 

10.00 

B3 V 

0.40 

-2.31 

HD 137569 

7.86 

B5ID 

0.40 

-5.63 

HD 62542 

8.04 

B5 V 

0.33 

-4.42 

HD 224424 

8.10 

B1 lab 

0.94 

-2.36 


* Dispersion normalized 2175A feature strength compared to the average 
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ABSTRACT 


Low resolution IUE data obtained towards HD 62542, a B5 V star in a line of sight associated with a 
dark cloud in the Gum Nebula region near £ Pup (05 Iaf) and yVel (WC8 + 091), yields an extinction 
curve with an extremely broad and weak UV extinction bump and the highest far-UV extinction yet 
observed in the Milky Way, reaching a value of E(1200-V)/E(B-V)~15. The weak bump is very similar 
to die bump observed towards HD 29647, a highly reddened late B HgMn star in the Taurus dark cloud. 
Using HgMn comparison stars, we have derived a new curve for HD 29647 which we believe accurately 
represents the extinction in this line of sight The far-UV extinction toward HD 29647 is larger than 
normal but not as extreme as that observed for HD 62542. 

Using analytic fitting techniques, we have derived parameters which describe the central position and 
FWHM for the weak bumps observed towards HD 62542 and HD 29647. In both cases, the bump 
central positions are shifted shortward of 2175 A with A. o =2110 A for HD 62542 and X 0 =2128 A for 
HD 29647. For HD 62542 and HD 29647, we find the bump FWHM to be 1.29 pm 1 and 1.62 pm' 1 
respectively as compared with 0.992 pm' 1 for die average of the stars in Fitzpatrick and Massa (1986). 
These two lines of sight provide the most extreme examples of anomalous UV bump extinction yet 
observed and reveal that there are occasionally substantial changes in the postion and width of the bump 
which may be produced by environmental factors. 

The dense region that gives rise to the extinction observed towards HD 62542 is exposed to both 
intense UV radiation and stellar winds from £ Pup and the yVel system. Strong shocks are undoubtedly 
exerting influence on the grains in this region. In contrast, the region responsible for the observed 
extinction towards HD 29647 is representative of a typical cold and quiescent molecular dark cloud. 
There is no evidence of any shocks in this line of sight The similarity of the observed UV extinction 
bumps in these two regions coupled with their extreme positions and widths and highly contrasted 
environments provide important new hints to the origin of the material responsible for the 2175 A bump. 
The extreme nature of the far-UV extinction observed toward HD 62542 may also help to constrain the 
range of possible explanations of the far-UV rise in interstellar extinction. 
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L INTRODUCTION 

The ubiquitous 2175 A interstellar extinction feature is by far the strongest known interstellar 
spectroscopic feature. For lines of sight with normal interstellar extinction [e.g. Seaton (1979), Savage 
and Mathis (1979)] and E(B-V)=1, the 2175 A interstellar extinction feature or bump as it is frequently 
called has an equivalent width of about 1000 A. By comparison, the same line of sight will typically 
have a neutral hydrogen column density of about 5X10 21 atoms cm' 2 which will produce an HI Lyman a 
line having an equivalent width of about 50 A. Among other potential dust features, the strongest 
measurable diffuse interstellar bands are typically less than 2 A (Herbig 1975). The origin of the 2175 A 
feature must rank as one of the most important unsolved spectroscopic problems involving the 
interstellar medium. 

Recently, serious attempts have been made to measure reliably the properties of the 2175 A feature in 
diffuse and cloud lines of sight through a parameterization of the bump's strength, shape and position 
(Massa and Fitzpatrick 1986; and Fitzpatrick and Massa 1986 ). Fitzpatrick and Massa (1986) fit the 
bump profiles for extinction along 45 lines of sight including a number exhibiting peculiar extinction 
(e.g. Her 36 and 6 Qri). Their results indicate that the bump position is very stable with an extreme 
deviation from 2175 A of only ±17 A. However, the bump width (i.e. FWHM) shows large variations 
of anywhere from 0.77pm' 1 (360 A) to 1.25pm' 1 (600 A) and the bump strength normalized to E(B-V) 
also shows large variability. The bump width seems to correlate roughly with crude indicators of line of 
sight dust density as measured by E(B-V)/r, where r is the distance to the reddened star. As discussed 
by Savage (1975), the very restricted range in variation of the bump position places very severe 
constraints on the possible range of parameters of the particles responsible for the feature. For recent 
discussions of this problem see Hecht (1986) and Mathis (1987). 

While the centroid of the 2175 A absorption feature for dust in the diffuse interstellar medium seems 
remarkably constant, it is noteworthy that several objects with circumstellar dust exhibit broad ultraviolet 
absorption features at wavelengths significantly different than 2175 A. A 2400 A absorption feature is 
seen toward R CrB stars and the carbon rich planetary nebula Abell 30 (Hecht et al. 1984; Greenstein 
1981). Some of the B and early A stars with circumstellar dust studied by Sitko, Savage, and Meade 
(1981) have extinction in the middle ultraviolet for which the 2175 A bump is either weak or missing 
entirely and the underlying extinction often contains shoulders or slope changes in the region 2500 and 
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2000 A (see Hecht and Nuth (1982) for additional discussions of these objects). Unfortunately, a 
serious complication in interpreting many of the ultraviolet circumstellar dust measurements is the 
presence of absorption from the circumstellar gas which can be confused with absorption by the dust 

In order to probe the extreme range of fundamental properties of the 2175 A feature, we have initiated 
an observing program with the International Ultraviolet Explorer (IUE) satellite to obtain measures of 
ultraviolet extinction toward stars known through broad band UV photometry to have highly abnormal 
UV extinction. The object selection is based on the catalog of UV extinction measures provided by ANS 
satellite UV photometry (Savage et al. 1985). The results for one of the objects we observed 
(HD 62542) are so interesting, we are reporting on the data here separately. Furthermore, the nature of 
the extinction bump in this line of sight is very similar to the feature observed towards HD 29647 which 
was presented in preliminary form by Snow and Seab (1980). We are therefore including a 
re-evaluation of that curve here. Finally, we have undertaken an analysis of the properties of these two 
extremely peculiar extinction curves utilizing the parameterization scheme of Fitzpatrick and Massa 
(1986) and discuss the results in context with theirs and the environments in which these interesting lines 
of sight reside. 


EL OBSERVATIONS AND REDUCTIONS 

The observed objects are listed in Table 1. The UV data for HD 62542 were obtained by us. The 
remaining UV data were obtained from the IUE archives. The sources for the photometry and spectral 
classification are also listed. All of the UV data were reduced at the Midwest Astronomical Data 
Reduction and Analysis Facility (MADRAF) at the University of Wisconsin, Madison. With the 
exception of HD 62542, all of the long- wavelength data were obtained with the LWR camera. In order 
to standardize the observations, all LWR images were corrected for LWR sensitivity degradation using 
the results of (Havel et aL (1986). No such correction is currently available for the LWP camera. Also, 
little information is available on the uncertainty associated with comparing LWP and LWR images. 
Clayton and Fitzpatrick (1987) have compared observations of a few stars obtained with both the LWP 
and LWR cameras and find no significant systematic differences between the two cameras. However, 
Harris and Cassatella (1985) have shown that apparent problems in the intensity transfer function (TTF) 
seem to result in systematically higher fluxes in the LWP camera at the overlap with the SWP camera for 
weak exposures. The shift in the spectrum appears to be confined to the region 1900 - 2100 A. For 
reasonably exposed spectra (data number=DN £ 130) the shift appears to be less than 10%. For our 
three LWP images of HD 62542, the exposure at 2100 A ranged from 130 to 1 80 DN (camera saturation 
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occurs at 255 DN). We feel that any systematic errors in our combined LWP data are less than 10% (see 
figure 1) and probably only affected the region below 2100 A. Recent tests comparing a new LWP ITF 
(Cassatella and Lloyd 1987) with the one used with our data indicate that a decrease in the flux below 
2200 A of order 3-5% relative to the region X>2300 A occurs if the new ITF is used. We shall discuss 
the implications of this in § IVa. 


a) HD 62542 

The reduced flux data for HD 62542 is shown in Figure 1. Also shown are the data for the two B5 V 
comparison stars, 23 Cyg and p Aur, listed in Table 1. The spectrum of HD 62542 was constructed by 
splicing together the individual observations listed in Table 1. These observations included both normal 
and heavily exposed images which were necessary in order to obtain improved signal-to-noise in the 
region of the bump as well as for X< 1600 A. For the combined spectra for each camera, the saturated 
portions were removed and the individual images were binned together. Finally, the combined SWP and 
LWP images were spliced together in the overlap region 1900 A - 2000 A. The only apparent 
discrepancy in the fit occurs in the region 1900 A - 1950 A where the combined LWP data appear to be 
about 5-10% higher than the SWP data. We subsequently gave the SWP data higher weight in this 
region when we combined the two images. 

The V and B-V photometry for HD 62542 is from Cousins et al. (1961) and Cousins and Stoy 
(1962). The B5 V spectral classification comes from Feast et al. (1955), Cousins and Stoy (1962) and 
Cousins et al. (1962). No peculiarities are noted in any of the sources which would lead one to 
question the spectral class assignment Unfortunately, no U-B color is available. However, Cousins et 
al (1961) did determine a (U-B) c color (measured on the Cape Observatory system). Although the two 
systems are generally not transformable, a crude transformation was constructed using a large number of 
reddened and unredded B0-A0 stars found in the references quoted above. The quoted (U-B) c color 
corresponded to a range of reddened U-B colors (using B-V=0.17) appropriate for spectral types B4 V 
to B6 V. 

The UV spectrum of HD 62542 was placed on a grid of UV comparison stars obtained from the IUE 
Atlas of Wu et al. (1983) in order to obtain the best UV spectral match. The grid ranged from B1 V to 
B9 V and included type in stars as well. The stellar features used for comparison included Si II 
(X1260), Si HI (XI 300), Si II (X1306), C U (X1334), Si IV (X1400), and the luminosity-sensitive Fe m 
features at X1900. For a complete discussion on the behavior of these features with temperature and 
gravity, see Panek and Savage (1976) and Koomeef and Code (1981). Types earlier than B3 were 
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eliminated due to moderate to strong (£ 10 %) mismatch for nearly all of the spectral lines, especially Si 
II (XI 260) and Si HI (X1300). For type B6, the Si II lines (XI 260, X1306) are noticeably stronger and 
grow rapidly for later types (40 % stronger than the lines seen in HD 62542 by type B7). For types 
B3-B5, no significant differences in the spectral line strengths longward of La are apparent at the 
resolution of our data. There is some stellar absorption in the region 1 190-1200 A that rapidly increases 
in strength for types later than B3. From this, type B5 seems to yield the most suitable comparison. 
However, the data in this region are weak and thus not extremely reliable and so we cannot totally 
discount types B3-B4 V on the basis of the UV observations. 

In establishing a UV luminosity classification of HD 62542, we note that the Fe HI features at 1900 A 
are very luminosity sensitive and were sufficient to confine the luminosity class to the range III - V. For 
types B3 - B5, the differences in the feature strengths between class HI and V objects are relatively small 
and so we cannot be definitive about the luminosity class from the IUE data alone. However, based on 
the fact that all three of the optical classification spectrograms indicate class V, we feel that it is justifiable 
to exclude giants from the discussion. 


b) HD 29647 

Figure 2 displays the reduced flux data for HD 29647. The optical photometry in Table 1 was 
obtained from Crutcher (1985) and references therein. As with HD 62542, the data shown in Figure 2 
are the result of combining the SWP and LWR images listed in Table 1. LWR 5641 was saturated in 
the region X>2450 A in order to produce reasonable S/N around the bump (Snow and Seab 1980) and 
so the weaker spectrum LWR 8532 was fitted in the region 2400 A - 2450 A. 

Straizys, Wisniewski, and Lebofsky (1982) and Straizys, Cemis, and Hayes (1985) report that 
HD 29647 shows abnormally strong Hg and Mn lines placing it in the peculiar subclass of HgMn stars 
(for a good discussion see Wolff (1983) and references therein). Straizys, Cemis, and Hayes (1985) 
also give log g = 3.7 indicative of class IV. Close examination of optical spectrophotometry obtained 
by one of us (JAC) indicates that the Balmer jump is too large for a main sequence star of the same 
temperature as inferred from the strength of the hydrogen Balmer lines. We estimate that the gravity is 
probably indicative of a late B giant. This is supported by Crutcher (1985) who indicates a luminosity 
ofm-IV. 

The optical line-blanketing of HgMn stars is small (Wolff 1967), and consequently the optical 
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continua of HgMn stars do not appear significantly different from normal stars of the same temperature. 
Hence, the U-B and B-V colors of HgMn stars do not appear to deviate from their normal star 
counterparts (Cameron 1966). However, there is evidence that the UV continua of HgMn stars deviate 
from their normal counterparts. Using OAO-2 data, Leckrone (1973) analyzed a number of 
abundance-peculiar A and B stars in the UV, including a number of HgMn stars. He found that HgMn 
stars tend to be underluminous in the UV compared to normal stars of similar temperature. Depending 
on the temperature, the differences range from 0.3 to 0.5 mag in the region 1400 - 2400 A and are be 
wavelength dependant. 

The treatment of HD 29647 and its subsequent peculiar extinction in the literature has generally been 
handled from the perspective of a normal late B main sequence star. The primary extinction curve for 
HD 29647 derived by Snow and Seab (1980) used a class la comparison star which we consider to be 
inappropriate. The data we present provide complete coverage in the UV region (the Snow and Seab 
data had a gap in the region 2400-3100 A due to IUE detector saturation). Finally, the type V 
comparison stars used by Snow and Seab (1980) were obtained from the TD1 data of Nandy et al. 
(1976). The known differences between the TD1 and IUE absolute calibrations (Bohlin et al. 1980) were 
apparently not allowed for in calculating the observed extinction curves. Also, based on the above 
luminosity classification discussions, comparison of HD 29647 to normal class V stars is inappropriate. 
We feel that the best comparison stars for HD 29647 should come from the HgMn subclass. Only a few 
such stars have been observed with IUE. Two such stars, k Cnc and x> Her, are listed in Table 1 and 
appear in Figure 2. Placement of the UBV data of these two stars on a two color diagram for class in 
stars (Fitzgerald 1970) indicates that the U-B and B-V colors are essentially the same as those for 
unreddened normal stars. The spectrum of HD 29647 seems to indicate additional line absorption not 
seen in either HgMn comparison star. The most notable differences occur at about XI 470, X2000, 
X2140, and X2500 A. The feature near X2500 A may be produced by strong Mn II resonance lines. 
Even with these absoiption line differences, we feel that from the limited UV work on these stars, the 
intrinsic continua is best matched with HgMn stars of the appropriate temperature range. 


IIL EXTINCTION CURVES: DERIVATION AND UNCERTAINTY 

The extinction curves presented here were derived by the standard pair method. Briefly, this involves 
comparing the colors m(X - V) of the reddened program star with the dereddened colors of a lightly 
reddened or unreddened comparison standard. The derived curves are then normalized by the color 
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excess of the reddened program star. The comparison stars for HD 62542 were taken from Wu et al. 
(1983). For HD 29647, two late B HgMn stars were chosen from the IUE archives. Below, the 
derived extinction curves and possible uncertainties are discussed separately for each star. 

a) HD 62542 

The extinction curve derived for HD 62542 using a B5 V comparison star appears in Figure 3. Also 
plotted is the average Milky Way extinction curve of Seaton (1979). The specific curve shown for HD 
62542 was derived using p Aur because the intrinsic reddening for this standard is the smallest of the 
two comparison stars (E(B-V)=0.01). However, the curve using 23 Cyg (E(B-V)=0.03) is essentially 
the same (see below). The gap in the extinction curve at X' 1 ® 8.2 pm" 1 corresponds to the position of 
La which was strongly contaminated due to the effects of geocoronal La in the observation of HD 
62542. 

Considering the relatively small E(B-V), the curve in Figure 3 is remarkably smooth. The maximum 
scatter in the normalized curve for X' 1 > 5 pm" 1 is of the order of ±0.25 mag which corresponds to a 
relative uncertainty in the flux ratio of less than 10%. No noticeable effects of absorption line mismatch 
appear above the noise, which implies the lines have been matched to better than 10%, a fact which 
confirms that our choice of comparison star is probably reliable. For X" 1 < 5 pm" 1 , the scatter is 
somewhat larger being about ±0.4 mag in the region where the LWP camera has a low sensitivity and 
about ±0.3 mag elsewhere. 

As pointed out in § n, spectral types B3-4 V cannot necessarily be ruled out on the basis of the UV 
spectral classification alone. The relative scatter discussed above also applies to extinction curves 
determined with B3 V comparison stars despite the fact that the normalized extinction curves are slightly 
different However, we suspect that type B3 V might be inappropriate for the following reasons. For 
types B2 V - B5 V, the B aimer jump rapidly increases at a rate of nearly 10% per spectral class (Kurucz 
1979; Jacoby et al. 1984). Use of a B3 V comparison star results in a normalized curve that is about 
+0.4 mag above the average curve at X _1 =3.2 pm" 1 and about +0.5 mag above the curve for HD 62542 
shown in Figure 3. Although there is no reason to expect our normalized curve to match the average 
curve at the overlap with the optical, most extinction curves including those for peculiar lines of sight 
generally normalize at or below the average curve. Comparatively few properly normalized curves are 
found above the average curve at X'^3.2 pm" 1 and to our knowledge none are more than 0.3 mag. 
Although eliminating type B3 V on this account alone is admittedly uncertain, the effects we see are 
completely in keeping with a mismatch at the Balmer jump. 
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In order to fully understand the accuracy of the curve shown in Figure 3, one must consider all 
possible sources of uncertainty, both random and systematic. A complete discussion of the sources and 
magnitude of these errors are presented in detail in Massa, Savage and Fitzpatrick (1983) and Massa and 
Fitzpatrick (1986). Random errors arise in general from uncertainties in the visual photometry and 
zero-point photometric errors inherent in calibrated IUE data. For the normalized curve shown in 
Figure 3, use of the propagation-of-eirors analysis of Massa, Savage, and Fitzpatrick (1983) indicates 
that the random errors are of the order of ±0.2 mag at A,' 1 =3.2 pm' 1 , ±0.4 mag at A' 1 =6.5 pm' 1 and 
up to ±0.7 mag at A' 1 =8 pm' 1 . Across the entire UV curve, this corresponds to a relative uncertainty of 
only ±0.5 mag. [We have assumed a(B-V)=0.015, a(V)=0.015, and the systematic shifts appropriate 
for the LWP camera to be similar to those discussed for the LWR camera (see Bohlin et al. 1980).] 

The major sources of possible systematic errors in our derived curve include comparison star 
temperature mismatch, improper dereddening of the comparison star, and errors introduced through 
possible ITF errors in the LWP camera data. As discussed in § II, this latter source will be confined 
primarily to the region around 2000-2100 A. For the moment we will ignore this effect. 

Using the convention of Massa, Savage, and Fitzpatrick (1983), the errors introduced through 
improper dereddening of the comparison star have the form [E(B-V) c /E(B-V)]x[k(A-V) c -k(A-V)] 
where E(B-V) C and E(B-V) are the color excesses of the comparison and program stars respectively, 
k(A-V) c is the normalized extinction curve appropriate for the comparison star, and k(A-V) is the 
assumed average extinction curve (e.g. Seaton 1979). For E(B-V) C =0.04 and Ak(A-V)=l, the error 
would be of the order of only 0.13 mag. Comparison of the curves derived by using 23 Cyg 
(E(B-V) C =0.03) and p Aur (E(B-V) C =0.01) shows that the two curves differ by less than ±0. 15 over the 
entire UV wavelength range. Given the low reddening of our comparison stars, we feel that the effect of 
this type of error in our curve will be less than 0. 1 mag and is consequently of little importance. 

Although we feel that B5 V is the best comparison star for HD 62542, as discussed in § II we are 
unable to completely discount type B3-4 V. However, for early BV stars, the effect on a normalized 
extinction curve of mismatching by 1-2 spectral classes is generally small. This results from the fact that 
the analytic expression for the mismatch error is somewhat self-cancelling in nature because it is 
proportional to the ratio of the terms Am(A-V) and E(B-V). Thus use of a hotter comparison star results 
in an increase in both Am(A-V) and E(B-V). For the extinction curve for HD 62542, we have examined 
mismatch errors empirically by comparing curves derived with both a B3 V and B4 V with the one 
shown in Figure 3 (B5 V). For the B3 V/B5 V case, Ak(A-V)=+0.5 mag over the region 3.2 pm' 1 < A' 1 
< 6.2 pm' 1 . Beyond A' 1 =6.2 pm' 1 , Ak(A-V) decreases monotonically to a value of +0.25 mag at A' 1 =8 
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pm' 1 . For the B4 V/B5 V case, Ak(X-V) increases from +0.3 mag to +0.4 mag for 3.2 pm' 1 < X* 1 < 
5.5 jim' 1 and decreases again monotonically to +0.3 mag at X _1 =8 pm' 1 . For both cases, Ak(X-V) 
shows negligible wavelength dependance for X' 1 < 6 pm' 1 . As discussed in the beginning of § HI, these 
apparent zero-point shifts are totally in keeping with mismatch at the Balmer jump. Across the entire 
curve, the relative error due to mismatch in the worst case discussed here is less than 0.3 mag. 

Collecting the errors discussed above and noting that their contributions are independent of each 
other, we find that the mean error relative to X' 1 =3.2 pm' 1 for the curve shown in Figure 3 is of the 
order of ±0.3 mag. 


b) HD 29647 

The extinction curve derived for HD 29647 appears in Figure 4. Also plotted is the average curve of 
Seaton (1979). The comparison star used was x> Her although the curve derived using k Cnc is 
essentially identical. The data for X' 1 £8 pm' 1 are very noisy and have been eliminated. Overall, the 
scatter in the curve is quite small with the exception of those features described in § II that showed no 
counterpart in either comparison star. Compared to the limited curve derived by Snow and Seab (1980), 
the scatter seen in Figure 4 is noticeably less. We feel that this is in part due to the better spectral match 
obtained by using HgMn comparison stars. 

As can be seen in Figure 4, the normalized curve has a value that is +0.3 mag above the average 
curve at X' 1 =3.2 pm* 1 . This shift is not due to mismatch at the Balmer jump, but represents the true 
position of the curve. This is confirmed by examining the broad-band optical extinction curve of 
HD 29647 determined by Straizys, Wisniewski, and Lebofsky (1982) which has a value that is +0.3 
mag above the average curve at X' 1 =2.9 pm* 1 . An optical extinction curve derived from 10 A 
resolution spectrophotometry obtained by one of us (JAC) yields a curve that is +0.25 mag above the 
average curve at X* 1 * 3.1 pm* 1 . 

As with the curve for HD 62542, we estimate random errors using the propogation-of-errors analysis 
of Massa, Savage, and Fitzpatrick (1983). We find that these errors range from ±0.06 mag at 
X' 1 =3.2 pm* 1 to ±0.2 mag at X* 1 =8 pm* 1 corresponding to a relative uncertainty of < ±0.15 mag. 
Considering the size of the color excess towards HD 29647, E(B-V)=1.01, it is not surpising that the 
errors are this small. A similar situation also applies to the errors corresponding to improper 
dereddening of the comparison star and spectral mismatch. In the case of the former, we have assumed 
E(B-V)=0 for both HgMn comparison stars in our derivation of the extinction curve based on the 
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discussion of the colors of HgMn stars in § lib and so no reddening correction was applied to the 
comparison stars. If the colors for these stars are in fact peculiar and E(B-V)*0, then from the two color 
diagram we require that the amount by which both U-B and B-V deviate from normal colors be exactly 
balanced by the effects of reddening. However, we find such a situation to be unlikely. Thus, if the 
colors are indeed normal and there is any reddening, it is probably hidden in the uncertainty of the optical 
photometry and is consequently of the order of E(B-V)=0.015. In this context, an upper limit to this 
error can be found by considering the effect of not correcting for reddening at all. From this it follows 
that the relative error across the UV curve will be at most 0.06 mag. Even if E(B-V)=0.03, the effect is 
still rather small and consequently of little relative importance. 

Because of its peculiar nature, assessment of spectral mismatch errors for HD 29647 is somewhat 
more difficult than for normal stars. Examination of the data of Leckrone (1973) seems to indicate that 
HgMn stars show no more variance in their UV flux distributions from object to object than do normal 
stars of the same temperature. However, these data represent only a limited sampling of wavelengths in 
the UV and unfortunately, the two HgMn stars listed in Table 1 represent the only existing high quality 
low resolution IUE data of this subclass of objects with die appropriate temperature. As stated above, 
both comparison stars yield essentially the same curve within the random errors. We therefore have 
taken the approach of examining variances incurred if normal stars of similar temperature are used. Since 
we know that use of normal type V stars is certainly incorrect, such a procedure establishes a suitable 
worst case error estimate. For this purpose we have chosen to compare the differences, Ak(X-V), 
between the curve in Figure 4 and those derived by using a B7 V, B8 V, and B8 IB comparison stars. 
For both the B7 V and B8 V curves, Ak(X-V) « +0.2 mag for 3.2 pm' 1 < A,' 1 < 5 pm' 1 . For X' 1 > 5 
pm' 1 , Ak(X-V) rises steadily to a value of +0.6 mag at X' 1 = 8 pm' 1 . Because HD 29647 is luminosity 
class IV-IH, this initial shift of +0.2 mag at the overlap with the optical is in keeping with a mismatch at 
the Balmer jump. This effect is clearly seen in the optical extinction curve when type B7 V or B8V 
comparison stars are used. For the comparison with a B8 HI, Ak(X-V) » 0 for 3.2 pm' 1 < X* 1 <4 
pm* 1 , which supports the Balmer jump mismatch for the type V curves discussed above. For X' 1 > 4 
pm* 1 , Ak(X-V) decreases rapidly to -0.1 mag and stays constant till around X' 1 = 8 pm" 1 where 
Ak(X-V) * -0.15. 


From the above discussion, we estimate the mean errors for the curve shown in Figure 4 are < 
±0.2 mag and are probably dominated by die random errors. It appears that using a normal comparison 
star of type B8 III does not significantly affect the shape of the curve. It would seen that the most 
serious source of mismatch arises from using stars of the wrong luminosity class. We feel that the curve 
in Figure 4 is well determined and thus represents an accurate assessment of the extinction towards this 
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line of sight. 


IV. EXTINCTION CURVES: PARAMETERIZATION 

From examination of Figures 3 and 4 it is easy to see that the normalized extinction towards these two 
lines of sight is quite peculiar, particularly in the region of 2175 A. In order to understand the properties 
of the dust in these two regions, it is helpful to analytically quantify the properties of their extinction 
curves particularly in the region of the 2175 A feature. To this end, we have chosen to use the analytic 
parameterization scheme developed by Fitzpatrick and Massa (1986). This procedure involves 
performing a least squares fit to the observed extinction using a linear (in A.* 1 ) background, a bump 
'profile' function over the range 3.2 pm' 1 £ A* 1 <, 5.9 pm' 1 , and a cubic function for A.' 1 > 5.9 pm' 1 . 

Savage (1975) showed, using OAO-2 data, that a linear background plus Lorentzian profile yields a 
good fit to the observed shape of the extinction bump. Seaton (1979) verfied this result using 
Copernicus and TD1 satellite data. More recendy, Massa and Fitzpatrick (1986) showed, using IUE 
observations of 33 stars in 5 open clusters, that such a procedure accurately codifies the shape of the 
extinction bump. Fitzpatrick and Massa (1986) have since performed parameterized fits for extinction 
observed in 45 lines of sight including a number exhibiting some of the most peculiar extinction 
observed in the Galaxy. However, they adopt a bump profile referred to as a Drude function which can 
be understood as an absorption cross-section under the assumption of the Drude theory of metals (see 
Bohren and Huffman 1983). They find that in addition to being easier to physically interpret than a 
Lorentzian, the Drude function also yields slighdy better fits. We have therefore adopted this approach, 
both for its validity and the fact that the results of Fitzpatrick and Massa (1986) supply a large data base 
to which are results can be reliably compared. The parameters determined from the bump fit include 
A^ -1 and y, the bump central position and width (i.e. FWHM) respectively and are discussed below. 
The function used to fit the observed extinction is given in Table 2. 

a) HD 62542 

The extinction curve for HD 62542 clearly represents one of the most extreme examples of extinction 
yet observed in the Galaxy. In the region of the bump, the extinction appears relatively broad and flat. 
In fact, if no prior knowledge of the shape of UV extinction is assumed, it is not clear that one would 
consider the region around 2200 A in Figure 3 to be a bump at all. In other words, the bump region is 
so broad and flat that one has to question if a bump is even present in the traditional sense. It is possible 
that what is observed is the underlying nature of UV extinction in the absence of the bump. 
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One possible reason why the bump region appears as a broad plateau is the rather large amount of 
far-UV (FUV) extinction present in this line of sight. Although not exceedingly steep, the FUV reaches 
values that are larger than for any curve so far measured in the Milky Way. Compared to the Milky Way 
average of 6.4 mag, E(1300-V)/E(B-V)=12 mag for HD 62542. By comparison, the value is 10 for the 
LMC 30 Dor region (Fitzpatrick 1986) while the value for the SMC is 14 (Prevot et al. 1984). For 
HD 62542, the extinction value at V 1 = 8.5 pm* 1 implies that the normalized extinction may reach 15 
mag at the shortest wavelengths accessible to IUE. 

Except for slope changes in the middle UV, which could be interpreted as an extinction bump, the 
observed extinction curve in Figure 3 is featureless. We will assume for the moment that a bump is 
indeed present in the general region around 2200 A. The results of the least-squares parameterization are 
shown in Figures 3 and 5a. In Figure 3, the fit to the extinction curve is overplotted on the observed 
curve. Also plotted are the individual components to the fit. As can be seen, the quality of the overall 
fit is good. In Figure 5a, the observed bump profile, found by subtracting the linear background and 
cubic FUV fit, is shown normalized to a peak value of 1 along with the determined Drude profile fit As 
can be seen, the observed bump is fitted quite well. Also plotted for reference is the equivalent Drude 
function appropriate for the average interstellar curve of Seaton (1979). 

Table 2 lists the results of the bump fit discussed above. The bump area, A^nAj/2y where A 3 is the 
bump normalization coefficient (see Table 2), is representative of the relative bump strength and is 
included for reference. Also listed are the extreme values of Xq** 1 and y found by Fitzpatrick and 
Massa (1986) along with the average for the 45 lines of sight they studied. Fitzpatrick and Massa (1986) 
empirically examined errors for the various bump parameters per unit E(B-V) that are expected for 
individual extinction curves. Because these errors were determined empirically from IUE data, they 
include the effects of random errors discussed in § IL Given that our fit shows no noticeable systematic 
errors, we have adopted their results as first order uncertainties. From their Table 2 we find 
CT(X o ~ 1 )=d0.02 pm* 1 , o(y)=±0.06 pm* 1 , and c(A)=±0.45 pm* 1 . The variance listed in our Table 2 
were derived by determining the bump parameters using both B5 V comparison stars. As can be seen, 
these empirical uncertainties are essentially the same as those derived from Fitzpatrick and Massa (1986). 
We also examined the bump parameters that result if spectral types B3 V and B4 V are used to derive the 
extinction curve. In both cases, A. 0 -1 and A were within the variances shown in Table 2. y was 
approximately 10% larger than the mean value listed in the Table. The only other potential effect on the 
bump parameters is the possible ITF problems in the LWP camera. As discussed in § II, the most 
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serious effect is a systematically higher flux of order 5-10% in the region 2000-2100 A. Decreasing the 
observed stellar flux by 10% in this region would increase the observed extinction by about 0.3 mag 
over the interval 4.76 pm' 1 < A* 1 < 5.00 pm' 1 . Such an increase would serve to shift the bump central 
position slightly to even shorter wavelengths but would not seriously affect the other parameters. 

Compared to the results of Fitzpatrick and Massa (1986), X 0 _1 and y for the bump observed towards 
HD 62542 represent the most extreme values yet observed. In terms of wavelength, X o =2110 A 
compared with ^=2158 A (0 Ori), the most extreme value found by Fitzpatrick and Massa (1986). 
Because the bump that we observe has the same functional shape as those seen in other lines of sight, it 
is tempting to assume that it is caused by the same material. However, because its central position is so 
removed from the normal position of 2175 ±17 A, we must consider that it could arise from some 
underlying component that is normally hidden by the standard bump. Considering the general invariance 
of observed bumps, the existence of such an extreme case makes this a reasonable possiblity to explore. 


b) HD 29647 

Since its discovery (Snow and Seab 1980), no extinction bump has been observed that is as weak as 
the one seen towards HD 29647. In fact, the extinction curve is often refered as being without a bump, 
although Figure 4 clearly indicates that a bump-like structure is present. In some respects, the curve in 
Figure 4 is quite similar to the one for HD 62542 shown in Figure 3. Both curves are characterized by 
weak bumps and FUV extinction that is higher than average. This general similarity is one of the major 
motivating factors behind our re-analysis of the extinction towards HD 29647. 

The results of the parameterization of the extinction curve for HD 29647 are shown in Figures 4 and 
5b. As with the HD 62542 curve, the computed extinction curve fit is overplotted on the observed 
curve. Also shown are the three components to the fit. It should be noted that the portion of the curve 
for X' 1 > 6 pm' 1 is not fitted well. There appears to me more structure than can be accounted for by the 
functional form shown in Table 2. Relative to the fit shown in Figure 4, the structure is in the form of 
=0.3 mag of excess extinction for 6 pm' 1 <X* 1 < 7 pm" 1 and =0.3 mag less extinction for 7 pm* 1 
7.5 pm* 1 . In Figure 5b, the observed bump is plotted along with the best-fit Drude profile and the 
Drude profile appropriate for the average curve of Seaton (1979). As with the bump for HD 62542, the 
bump observed towards HD 29647 is fitted well by the Drude function. 

Table 2 lists the relevant bump parameters determined from the fit procedure. From Fitzpatrick and 
Massa (1986), we estimate the errors in the fit to be a(X o ' 1 )=±0.01 pm' 1 , c(y)=±0.02 pm' 1 , and 
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o(A)=±0.15 pm" 1 . The variances determined from using both HgMn stars appear in Table 2. With the 
exception of Xq* 1 , our empirically determined variances are considerably larger than those deduced from 
Fitzpatrick and Massa (1986). This is not surprising, however, considering how weak the observed 
bump actually is (see Figure 4). In analogy to weak and broad spectral features, centroids are generally 
much easier to determine than the FWHM or the equivalent width. In such a case, slight variations in the 
placement of the continuum can have large effects on the equivalent width because of significant 
information in the wings of the feature. 

Despite the larger errors, we consider the bump parameters for HD 29647 to be relatively well 
determined. We of course make this statement with some caution keeping in mind the peculiar nature of 
HD 29647. Because little information still exists on the UV behavior of HgMn stars, we cannot be 
totally sure that use of the comparison stars in Table 1 has not introduced systematic errors into our 
bump parameters. However, parameterization of the bump produced by using a normal B8 HI yielded 
results consistent with those listed in Table 2. The reddening to this line of sight is so large that the 
effects produced by extinction appear to dominate over any uncertainties associated with the stars 
peculiar status. We therefore accept the values listed in Table 2 as representative of the true extinction 
bump observed towards HD 29647. 

Like HD 62542, the bump parameters for HD 29647 appear to deviate significantly from any bump 
yet measured. Although the bump is not as shifted as the one observed towards HD 62542, X Q for the 
HD 29647 bump is still shifted shortward of the standard bump position by nearly 50 A. More 
impressive is the FWHM which, in units of pm" 1 , is 30% larger than the most extreme example (C Oph) 
found by Fitzpatrick and Massa (1986) and nearly 70% larger than the average. These values are so 
different from what is normally observed that it is again tempting to attribute the bump extinction 
towards HD 29647 as being due to a possible sub-component which is only observable in the absence of 
the normal bump component Clearly, such a situation implies special physical conditions which are 
generally not applicable to most observed lines of sight. In the following section, we explore this 
possibility for both the HD 62542 and HD 29647 lines of sight by examining the environments in which 
the peculiar dust extinction arises. 


V. INTERSTELLAR ENVIRONMENTS OF HD 62542 AND HD 29647 

The peculiar extinction bumps presented for HD 62542 and HD 29647 represent either a previously 
unobserved sub-component, or simply extreme modification of the normal bump carrier. The derived 
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parameters are so different than all previous observed bumps (including peculiar bumps) that the 
physical conditions in the regions in which the dust exists must be extremely rare. In this section, we 
examine the regions in which these lines of sight occur in order to explore any physical parameters that 
might be responsible for the observed extinction. 

a) HD 62542 

The line of sight to HD 62542 is in a region of the sky which places it in the direction of the Gum 
Nebula complex. In Ha emission, the complex appears roughly circular with a diameter of about 36° 
(Sivan 1974) and is comprised of patches and filaments of emission consistent with the image of an 
expanding bubble-like nebula (Reynolds 1976a). The distance to the center of the nebula is estimated to 
be of the order of 400 pc (Brandt et al. 1971; Gott and Ostriker 1971) which translates to a diameter of 
roughly 250 pc. The origin of the complex remains a mystery although a number of scenarios have been 
proposed. These include a 'fossil Stromgren sphere' produced by the Vela supernova remnant 
progenitor (Brandt et al. 1971), a normal H II region being currently ionized by imbedded hot stars 
(Beuermann 1973; Wallerstein et al. 1980 ), a 10^ year old supernova remnant currently being heated by 
hot stars (Reynolds 1976a,b), and an expanding bubble created by ionization and winds from imbedded 
OB stars (Weaver et al. 1977). 

Regardless of it origins, die observed Ha emission appears to be consistent with ionization from two 
major sources namely £ Pup (05 Iaf) and y Vel (WC8 + 091) [Reynolds 1976a,b]. It also appears that 
some energy may also be supplied by stellar winds in the form of shocked gas as indicated by measured 
values of [S m/Ha (Chanot and Sivan 1983). The likely sources of these winds are the WC8 star y 2 
Vel which has a mass loss rate of 7.4-1 1.6 x 10* 5 M^yr and a terminal wind velocity of V*, = 1600 
km/s (Abbott et al. 1986; van der Hucht, Cassinelli, and Williams 1986) and £ Pup which has a mass 
loss rate of 3-7 x 10 -6 M^/yr and a terminal wind velocity = 2700 km/s (Olson and Castor 1981). 

In projected angular distance, HD 62542 is located about 10° west-southwest of the projected center 
of the nebula defined by the Ha emission and resides coincident with a ridge of dark matter. Assuming 
a normal value of total-to-selective extinction and adopting My = -1.0 (Johnson and Iriarte 1958) we 
estimate the distance to the star to be of the order of 400 pc which places it inside the expanding shell at 
about the same distance as C, Pup and y Vel. With respect to the sources supplying the energy to the 
complex, HD 62542 is located 4.5° southwest of £ Pup and 7° northwest of y Vel which at 400 pc 
translates to tangential distances of 31 and 49 pc respectively. A photograph of the region where HD 
62542 resides is shown in Figure 6. The directions to £ Pup and y Vel are indicated in the Figure. As 
can be seen, the eastern edge of the ridge is extremely well defined giving the impression that matter has 
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been swept up from the direction of £ Pup and y Vel . Examination of the reddening of B stars located 
east of the ridge and to the southwest and northwest of C, Pup and y Vel respectively yields typical values 
of E(B-V)<0.05 suggesting that this is a cavity relatively void of dust and gas. This is confirmed by the 
Ha photograph of Sivan (1974) which shows this region void of emission. (TTiis cavity is rimmed by 
relatively bright emission but it is difficult to determine if we are actually looking at a vacant volume 
surrounded by Ha emission or simply a hole in the foreground portion of the expanding shell.) 

Superposition of the Ha photo on Figure 6 shows the ridge to be well defined by bright Ha emission 
with a relatively sharp cut-off at the eastern edge and more diffuse emission to the west. Faint emission 
is also visible in Figure 6. A major source of the photoionization is undoubtedly radiation from C, Pup 
and possibly the 09 1 companion of y 2 Vel [Hua et al. (1983) report that the ionizing UV flux from WC8 
stars would be negligible]. It is quite plausible that shocks produced by winds from y 2 Vel ( 1600 km/s) 
and C Pup (2700 km/s) are also producing some of the observed emission. The sharp eastern 
boundary, coupled with the ridge appearing bowed in regions of apparent lower density also seems 
consistent with this picture. Interaction of strong winds with matter in this region is supported by the 
presence of numerous cometary-shaped dark globules to the southwest and northeast of the ridge, the 
tails of which point in the general direction of the wind flow from y 2 Vel and £ Pup. The eastern-most 
portion of the ridge in which HD 62542 is located shows characteristics indicative of a cometary-like 
globule in formation. Figure 6 clearly shows this to be the densest portion of the ridge. To the north of 
HD 62542, the ridge appears to show the effects of turbulence which could be caused by shearing as the 
wind from y 2 Vel flows along the edge of the dense material. Finally, the inset in Figure 6 shows the 
location of IRAS point sources observed in the 12 and 60 p bands in the general direction of HD 62542. 
As can be seen, the majority of the point sources are coincident with the densest portion of the ridge. It 
is tempting to imagine that shocks might have triggered local collapse and possibly star formation. 
However, close examination of Figure 6 shows the matter in this part of the ridge to be rather clumpy 
and it is quite possible that what we are seeing are clumps of warm dust with peak emission in 60 p 
band. 

From the data of van der Hucht, Cassinelli, and Williams (1986), pressure balance between the wind 
from y 2 Vel and the ambient medium will be achieved when d 2 nT=7.8xl0 7 pc 2 cm* 3 K where d is the 
distance from y 2 Vel in parsecs, n is the ambient density in cm* 3 and T is the ambient temperature. For 
nominal diffuse interstellar conditions (T®80 K, n«20 cm* 3 ), balance will be achieved for d~360 pc. Of 
course, this estimate is an upper limit since it does not take into account the increase in density that 
would occur in front of the ’snowplow’ in the case of isothermal expansion. However, d will certainly 
be > 50 pc for the ambient conditions assumed above and so the wind could easily have evacuated the 
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region to the east of the ridge. For pressure balance to occur at a distance d=50 pc from t 2 Vel, nT« 
3xl0 4 cm’ 3 K. For a local gas temperature of 100 K, n«300 cm’ 3 behind the ridge is sufficient to 
balance the wind. From this, we can conclude that the ridge seen in Figure 6 could easily remain well 
inside the Ha shell at a distance of 50 pc from y 2 Vel. (We have used y 2 Vel as an example because the 
mass flux produced by the wind is about a factor 10 greater than that produced by £ Pup. Inclusion of 
the wind from £ Pup would not effectively alter the above discussion.) 

If the ridge consists of a fair fraction of swept-up material from the cavity, then shock processing and 
sputtering will have had a profound effect on the observed grains. Some fraction of the observed grains 
may even have come from dust formed in the wind near the star. Even if material in the cavity has added 
little new matter to the ridge, shock heating as well as cloud compression must have modified the grains 
as seen by the extinction shown in Figure 3. This processing may include destruction as well as 
coagulation in the denser environments. 


b) HD 29647 

HD 29647 is a highly obscured late-B star in Taurus. With E(B-V)=1.03, HD 29647 resides in a line 
of sight significantly more obscured than HD 62542 which has E(B-V)=0.33. Crutcher (1985) has 
provided a comprehensive discussion of the properties of the region in which HD 29647 resides. It is 
situated within the area of Heiles cloud 2 and very close to the Taurus Molecular Cloud 1 (TMC 1), 
which is a very dense molecular clump within cloud 2. Crutcher's distance estimate for HD 29647 of 
170+20 pc suggests that the star lies beyond the Taurus daik cloud complex which lies at 140±10 pc 
(Elias 1978). This view is supported by the appearance of the Palomar Sky Survey prints which show 
no indication of nebulosity which we would expect to see if HD 29647 were imbedded in the dense 
interstellar region. In addition, the K and L photometric data of Straizys, Wisniewski, and Lebofsky 
(1982) show no evidence of infrared excess which is occasionally seen towards stars imbedded in dense 
matter (e.g. p Oph). Finally, Crutcher (1985) finds that the column density derived from the optical CN 
absorption lines is in exact agreement with the column density derived from the millimeter CN emission 
lines which would only occur in the star was entirely behind the cloud. 

The Palomar Sky Survey prints of the region around HD 29647 show no evidence of the very 
pronounced stellar-interstellar interaction which appears to be affecting the material in the line of sight to 
HD 62542 (see § Va and Figure 6). Other indications of the quiescent nature of the HD 29647 line of 
sight come from the optical absorption line data of Crutcher (1985). Elitzur and Watson (1978, 1980) 
have proposed that CH + is formed in shock heated gas via endothermic reactions of C + with H 2 . Such 
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a model has proven to be necessary because CH + is expected to be rapidily destroyed through reactions 
with electrons and H 2 and in many cases, the predicted ratio of NCCHyNCCH 4 ") through equilibrium 
chemistry yields values as much as two orders of magnitude larger than the observations (Federman 
1982). For the line of sight to HD 29647, the data of Crutcher (1985) yields N(CH)/N(CH + ) >17. For 
HD 147889, a star located in the p Ophiuchi dark cloud with a reddening similar to HD 29647, 
N(CH)/N(CH + ) ~ 2.8 (Cardelli and Wallerstein 1986). The radial velocities of the CH and CH + 
towards HD 147889 differ by 1.2 km/s which also conforms to shock formation theory (these results 
are also indicated for the general Ophiuchus complex; see Meyers et al. 1985). No such differences in 
velocity are found for the lines observed towards HD 29647. This, along with the large ratio of 
N(CH)/N(CH + ) implies that shock production of CH + is negligible towards HD 29647. Thus, while 
shock processing may influence the extinction towards HD 62542 (unfortunately, no information is 
currendy available on the optical absorption), the process probably does not play a role for the extinction 
observed towards HD 29647. The line of sight to HD 29647 suggests that the extinction is primarily 
produced in a relatively quiescent but very dense interstellar region. The diagnostic radio data of 
Crutcher (1985) reveal a region with very strong molecular lines and imply a high density of 800 cm' 3 
and a low kinetic temperature of 10 K, values typical of dark clouds. 


VI. DISCUSSION 

Figure 7 compares the extinction toward HD 62542 and HD 29647 with the average extinction in the 
Milky Way (from Seaton 1979) , the Large Magellanic Cloud (LMC) and the 30 Doradus region of the 
LMC (both from Fitzpatrick 1986). The extinction curves for HD 62542 and HD 29647 are among the 
most anomalous so far observed for dust in diffuse or dense clouds. The sight lines to both objects 
exhibit extinction curves with middle ultraviolet extinction bumps which are very weak, very broad and 
shifted in wavelength from the normal position of 2175 A. In addition, both sight lines have larger than 
average FUV extinction with the curve for HD 62542 having the largest FUV extinction ever recorded 
for a Milky Way star. In fact this curve even exceeds the high FUV extinction seen toward the 30 Dor 
region in the LMC. The discussions of these two unusual extinction curves which follow will consider 
both the ultraviolet bumps and the FUV component of extinction. 
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a) The Extinction Bump 

Although the ultraviolet extinction bump is the strongest spectroscopic feature observed in the 
interstellar medium, its origin is not established. However, because of the great strength of the bump, it 
follows that the material responsible for the bump must be very abundant. A limit to the required 
abundance can be obtained by assuming the absorbing atoms (in either atomic, molecular or solid state 
form) have an absorption oscillator strength or f value of unity . The standard relation between 
absorption column density, N, equivalent width, , wavelength , X , and oscillator strength , f, for an 
absorption feature which is not saturated is, IX - (k ) NXf. In the case of the 2175 A 

feature associated with the average interstellar extinction, « 1000 A when E(B-V) = 1.0. The 
column density of absorbers having f = 1 required to produce such a large absorption feature is N(bump) 
** 2.4x10*7 atoms cm'^. For a normal gas to dust ratio (Bohlin, Savage and Drake 1978) , the total 
hydrogen colum density for E (B-V) = 1 is N(H) = 5.8x10^* atoms cm'^. Thus, the required 
abundance relative to H of the bump material for an oscillator strength of unity is 4.1x10'^. This 
number should be considered as the minimum abundance required since most atomic, molecular or solid 
state transitions have f < 1. Among the abundant heavy elements that might plausibly produce 
interstellar solids, only C, N and O are comfortably more abundant than this constraint Mg, Si and Fe 
have abundances roughly comparable to the constraint Thus, if some solid state form of Mg, Si or Fe is 
responsible for the bump, then essentially all of the interstellar Mg, Si, or Fe must be tied up in that solid 
and the f value of the absorption must approach unity. Although not stated in these terms, this 
abundance constraint has forced most theories of the 2175 A bump to concentrate on processes arising in 
solids composed of the most abundant elements. 

Observations providing new insights about the range of the spectroscopic behavior of the bump are 
potentially of great value to those experimental and theoretical studies which seek to determine the origin 
of the bump. The early measures of the extinction bump with the OAO-2 satellite (Savage 1975) revealed 
a feature with a profile showing a remarkable degree of constancy among a sample of 36 stars. The 
observation that the bump position and shape were constant or very nearly constant placed severe 
restrictions on the grain geometrical parameters if the feature was explained by classical scattering theory 
employing bulk optical constants. In fact the restrictions seemed so severe that Savage(1975) suggested 
that alternative explanations may be required unless the grains which exist in widely separated regions of 
space and under very different physical conditions have nearly identical size and shape distributions. 
These restrictions can now be relaxed slightly because the IUE measurements reported by Fitzpatrick and 
Massa (1986) have revealed a 25 % variation in the width of the ultraviolet extinction bump. However, 
for the Fitzpatrick and Massa sample of 45 extinction curves the bump exhibited an extreme position 
variation of only ±17 A or ± 0.8 %. For a summary of the Fitzpatrick and Massa results along with 
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extinction parameters measured for HD 62542 and HD 26947 see Table 2. Both HD 62542 and HD 
29647 exhibit bumps with position and width parameters which lie outside the extreme limits found by 
Fitzpatrick and Massa. 

General reviews of the physical nature of the types of solid state processes which might produce a 
bump like feature are found in Gilra (1972) and Savage (1975). A number of factors can influence the 
position and width of an extinction bump if it arises from classical scattering and absorption by 
interstellar grains . Some of these factors include: 1. particle optical properties which change as the 
physical and thermal environment of the grains change; 2. particle size; 3. particle shape; 4. presence or 
absence of coatings or the existence of grain surfaces which have optical properties which differ from 
those of the grain interior. The sensitivity of the bump position to these factors depends on which of the 
many possible physical processes actually produces die bump. 

In recent years, most of the theories for the production of the bump have focused on the process of 
plasma oscillations in small grains of carbon which are usually considered to be in the form of graphite. 
The sensitivity of the positions of extinction bumps to particle geometric properties is illustrated in Table 
3 where we list the values of the central position of extinction bumps for a variety of particle sizes and 
shapes computed with classical Mie scattering theory and the bulk optical constants of graphite from Taft 
and Phillips (1965). As can be seen, a rather modest modification of the grain geometrical 
characteristics produces large changes in the position of the computed bump. Table 3 also lists the 
expected positions of extinction bumps produced by amorphous carbon and and amorphous carbon 
annealed at 2200 C (see Hecht 1986). From these numbers we see that the degree of disorder in the 
solid state form of matter in the absorbing solid also plays an important role in det ermining die solid's 
optical constants and hence bump position. 

Various suggestions have been made to address the problem of the stability of the bump positiori. In 
the diffuse interstellar dust model of Mathis, Rumpl and Nordsieck (1977, hereafter MRN) which 
employs a power law size distribution, the bump is mosdy produced by plasma oscillations in graphite 
particles in or near the Rayleigh limit For such particles the position of the bump is relatively insensitive 
to changes in the particle sizes. According to MRN, the bulk optical properties of graphite as measured 
by different experimental groups exhibits large enough of a range so that spherical graphite particles in 
the small size limit may indeed have a bump centered near the average observed value of 2175 A . An 
alternate solution has been advanced by Hecht (1986) who has considered the implications of the 
Fitzpatrick and Massa (1986) measurements. In reviewing the many possiblilites, Hecht concludes that 
the 2175 A feature is caused by a population of small hydrogen free carbon grains. He notes that large 
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carbon grains (which will have lower temperatures than the small grains) will likely retain a surface layer 
of attached hydrogen atoms and that this coating will produce a substantial suppression of the bump. 
Hecht also proposes that the width of the bump .which changes considerably from one sight line to 
another, is probably influenced by temperature effects and size effects among the population of small 
(5 to 50 A) hydrogen free carbon grains which he believes causes the feature. 

Within the framework of the Hecht model for the ultraviolet extinction bump, the weak bumps for 
HD 62542 and HD 26947 might be explained by environmental effects. In the dense interstellar 
environments for these two sight lines perhaps even the very small carbon grains have surface layers of 
hydrogen and therefore weakened bump strengths. The great breadth of the bumps would then have to 
be explained through size distribution and temperature effects. The shifted bump positions would also 
have to be explained by such effects. 

While the discussion above seems reasonable, it is important to emphasize that a number of other 
effects may also be operating to modify the extinction bumps observed toward HD 62542 and HD 
29647. Compositional differences, grain coagulation, grain shattering or the growth of icy mantles are 
just a few of the many possibilities. The icy mantle suggestion was advanced by Snow and Seab (1980) 
to explain the anomalous ultraviolet extinction to HD 26947. A search for the 3.1 pm H 2 O ice feature 
by Whittet et al.(1981) toward HD 26947 was unsuccessful, while a more sensitive observation by 
Goebel (1983) revealed an IR absorption feature near 3.1 pm with a strength of about 10 %. Scattering 
and absorption calculations for graphite grains coated with crystalline H 2 O ice predict that the resultant 
extinction bumps are weakened and shifted to longer wavelengths (Whittet et al 1981). The weak bump 
for HD 26947, however, is shifted to shorter wavelengths. If the weak feature seen by Goebel 
actually is produced by H 2 O ice absorption, the ice might therefore exist on another component of the 
interstellar dust. Thin H 2 O ice coatings on small grains are expected to produce absorption edges near 
A - * = 6.37 pm (Field,Partridge and Sobel 1965; Whittet et al 1981; Massa, Savage and 
Fitzpatrick 1983) . Such an ice edge having a strength of about 0.3 mag could exist in the extinction 
curve for HD 26947 shown in Figure 4. However, we note that the visibility of the U V H 2 O ice edge 
will be sensitive to such variables as the distribution of sizes of the particles being coated as well as the 
distribution in sizes of the coatings themselves. Detailed model calculations will be required to evaluate 
the true implications of the undulations in the HD 29647 extinction curve for > 6 pm 

Although H 2 O ice coatings would be expected to shift the 2175 A feature to longer wavelengths, it is 
possible to have other kinds of coatings which shift the feature to shorter wavelengths. Gilra (1972) 
discusses the effect coatings have on the position of plasma resonances such as those produced by 
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graphite. In the small particle limit for spherical particles, the plasma resonance will occur between the 
wavelengths where ej(X.) = —2 and Ej ( A ) = -2 e c (A,) , where e^(A) is the real part of the dielectric 
constant for the spherical particle being coated and e c (A,) is the real part of the dielectric constant for the 
coating material . In the case of plasma resonances in graphite, the feature will be shifted to longer 
wavelengths for e c > 1 and to shorter wavelengths for e c < 1. While H 2 O ice has > 1 in the 
wavelength region of interest , it is interesting that NH 3 ice has e c > 1 for A, > 2000 A but near 1990 A 
the NH 3 absorption edge will cause to become < 1 for A. slightly less than 1990 A. The presence 
of this absorption edge near the bump position might actually create a broadening of the bump profile. 
Among other plausible coatings it is noteworthy that metals such as Fe have e c < 1 in the relevant 
wavelength range. 

One of the most intriguing aspects of the extinction bumps observed for HD 62542 and HD 29647 is 
that they more closely resemble the bumps expected for small spherical graphite particles than any of the 
bumps previously observed. Reference to Tables 2 and 3 shows that the 'normal' observed bump with 
Aq = 2175 A and y = 0.99 |Xnf ^ occurs at the wrong wavelength and with too narrow a width to be 
compatible with classical scattering by graphite having the bulk optical constants as measured by Taft 
and Phillips (1965). Of course, it is easy to imagine that graphite in the diffuse interstellar medium has 
somewhat different optical constants and that those constants produce a bump for particles in the 
Rayleigh limit near the observed wavelength with the appropriate width. However, it is an interesting 
concidence that the bumps toward HD 62542 and HD 29647 have positions and widths which are very 
compatible with the theoretical expectations for small bare spherical graphite particles in or near the 
Rayleigh limit Perhaps the bumps we are measuring toward these two stars are 'normal' for interstellar 
graphite grains and the so called 'normal' bumps for diffuse cloud matter are representative of 
interstellar graphite grains which have somehow become modified. 

Finally we cannot rule out the interesting possibility that the bumps measured toward HD 62542 and 
HD 29647 reveal the underlying graphite component of the interstellar dust and that the bump seen at 
2175 A along diffuse cloud lines of sight is produced by some other process. 

b) The Far Ultraviolet Extinction 

The extinction curves for HD 62542 and HD 29647 are very unusual in the FUV. In interpreting 
the significance of these curves, it is important to emphasize that the results given in Figures 3 and 4 are 
presented as extinction normalized to E(B-V) = 1.0. In such a normalized extinction curve, high FUV 
extinction might arise either from an enhancement of that small particle component responsible for the 
FUV extinction or by a reduction in the effectiveness of the particles responsible for extinction between 
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the B and V photoelectric bands. 

In the MRN interstellar dust model the FUV extinction is provided by intermediate graphite and small 
silicate grains. Thus an enhancement in the relative numbers of these grains compared to that component 
of the dust producing E(B-V) extinction could explain an enhanced FUV extinction. Seab and Shull 
(1983) showed that an attractive theory for providing such an enhancement is shock processing of 
interstellar grains. To test this theory, Seab and Shull modeled the shock processing of an MRN dust 
mixture and found significant processing for shock speeds as low as 40 km s'*. The calculation reveals 
that grain-grain collisions and non-thermal sputtering are the dominant destruction mechanisms and these 
processes destroy more large grains than small and destroy more silicates than graphite. The modeling 
reveals that the large silicate grains which mostly provide B-V reddening are the most susceptible to 
destruction. The result of the shock processing of a MRN mixture according to Seab and Shull is 
therefore a reduction in the value of E(B-V) for a given mass of dust which shows up as a strong 2175 A 
bump and large far-UV extinction in extinction curves normalized to E(B-V) = 1.0. Unfortunately , the 
basic input physics required for carrying out the theory of shock processing of interstellar grains appears 
uncertain. Strazzulla Barata and Magazzu (1985) have criticized the Seab and Snow shock processing 
theory. By examining the laboratory studies of low energy sputtering of solids, they conclude that 
graphite and silicate grains will be almost equally destroyed at shock velocities greater than 70 km s'* 
and at lower velocities graphite is preferentially erroded. Until the theoretical predictions are clearified, it 
appears dangerous to try to interpret extinction curve variations based on theories of shock processing. 

The FUV extinction toward HD 62542 is the highest so far found for a Milky Way star. The swept 
up nature of the region certainly suggests that shock processing may have played significant role in the 
early history of the dust which is now seen along the sight line to HD 62542. However, a number of 
other processes may also have occurred. A combination of shock processing followed up by the dust 
modifications which likely occur in dense interstellar environments as discussed in § Via will probably 
be required to provide a complete understanding of the unusual extinction toward HD 62542. 

The extinction curves for HD 62542 and HD 29647 go counter to the general trend of finding low 
FUV extinction curves for dust in dense interstellar environments (see Mathis 1987). In the case of HD 
62542, the swept up nature of the interstellar environment (see §Va) suggests that this may have 
occurred through shock processing of interstellar grains. However, in the case of HD 29647 such an 
explanation seems untenable since the local environment appears relatively undisturbed (see §Vb). In 
this latter case, a more plausible explanation is that grain growth through mandes or coagulation has 
reduced the (B-V) reddening efficiency of those grains which mosdy contribute to visual extinction but 
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has left the far ultraviolet extinction relatively unchanged with the result that the far ultraviolet extinction 
appears enhanced in an E(B-V) normalized extinction curve. The idea of icy coatings is consistent with 
the 3. 1 jim HjO ice feature observations of Goebel (1983) . Mantles which consist of ices of abundant 
species such as H 2 O, NH 3 or CH 4 can reveal their presence through their strong UV absorption edges 
near A,"* = 6.37, 5.03 and 7.41 [im respectively (Field, Partridge and Sobel 1965) and weak (Am = 
0.3 mag) edges near these positions could certainly exist in the observed curve for HD 26947. The 
sensitivity of the expected appearance of these absorption edges to assumed particle and coating 
parameters is demonstrated in the theoretical curves shown by Field et al . (1965) . 

The highly anomalous character of the UV extinction bump toward HD 62542 and HD 26947 
suggests another possibility for the unusual FUV extinction. Graphitic carbon exhibits a second plasma 
resonance which in the small particle limit peaks in absorption at about 800 A. Perhaps the physical 
effect which has altered the wavelength and breadth of the 2175 A feature has produced a major 
modification in the character of this second resonance such that the effects of this modification are visible 
in the FUV region accessible to the IUE. 

c) General Extinction Considerations for Stars in Dense Interstellar Environments 
The two UV extinction curves found for HD 62542 and HD 26947 should cause major concern 
about the validity of using an average extinction curve to correct UV energy distributions for the 
presence of interstellar dust. It has been common in recent years to assume that extinction curves with 
low FUV extinction (e.g. Orion like or p Oph like ) are generally applicable in dense environments. The 
curves for HD 62542 and HD 26947 reveal that there are cases where the extinction has the opposite 
behavior. The extreme nature of the curve for HD 62542 unfortunately reveals that the errors associated 
with using a standard curve can be enormous. For example, if a curve like that for 0 *C Ori (Bohlin 
and Savage 1981) is assumed to be valid but the actual curve is like that seen toward HD 62542 and if 
the standard process of "ironing out" the bump is followed then the corrected fluxes at 1250 A compared 
to the V band would be in error by an astonishing 1 1 mag for E(B-V) = 1.0. 

In addition to the standard problem of correcting astronomical observations for the presence of dust, 
there are important physical processes which are strongly affected by the dust properties. For example, 
interstellar chemical processes are influenced by the transfer of dissociating FUV radiation into dense 
interstellar regions. The exceedingly high FUV extinction found toward HD 62542 will likely play a 
major role in influencing the chemistry in the cloud toward this interesting star. Also, the presence of 
large numbers of small grains along this sight line may substantially modify the grain surface chemistry 
over that found in normal diffuse clouds. Follow-up high dispersion IUE and ground based 
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observations of narrow interstellar absorption lines in the spectrum of HD 62542 may therefore be 
exceedingly interesting. 

We acknowledge valuable discussions with J. Mathis regarding UV extinction and the theories of 
interstellar dust. D. Leckrone provided helpful suggestions regarding the UV energy distributions of B 
stars having peculiar abundances. Marilyn Meade provided assistance with the IUE data handling. We 
thank the IUE observatory staff for their help in acquiring and processing the satellite data. The final 
data reductions were processed at the Midwest Regional Astronomical Data Analysis Facility in 
Madison. This work was supported through NASA grants NAG- 186 and NAG-832. 
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TABLE 1 

REDDENED AND COMPARISON STARS 


Name 


V 

B-V 

SpType 

References 

IUE 

Images 

Aperture 

Exp.(sec) 

HD 62542 


8.04 

0.17 

B5V 

1 

SWP 28644 

L 

145 







SWP 29212 

L 

300 







SWP 29213 

L 

600 







LWP 8580 

L 

140 








S 

130 







LWP9077 

L 

60 







LWP 9078 

L 

240 

HD 34759 

p Aur 

5.22 

-0.15 

B5 V 

2.3,4 

SWP 15537 

L 

8.6 







LWR9868 

L 

7.7 

HD 188665 

23 Cyg 

5.14 

-0.13 

B5 V 

2,3,4 

SWP 15338 

L 

5.4 







SWP 15339 

L 

8.0 







LWR 11856 

L 

8.2 







LWR 12008 

L 

9.2 


HD 29647 


8.37 

0.90 

B8IV-mp 

1 

SWP 6557 

L 








LWR 5641 

L 








LWR 8532 

L 


HD 783 16 

k Cnc 

5.23 

- 0.11 

B7-8mp 

5,6,8 

SWP 15804 

L 

24 







SWP 15820 

L 

24 







LWR 12175 

L 

15 







LWR 12184 

L 

15 

HD 144206 

\) Her 

4.76 

- 0.11 

B9mp 

5,6,7 

SWP 4965 

L 

7 







LWR 4297 

L 

4 


(1) see text 

(2) Wu et aL (1983) standard star 

(3) Lesh (1968) 

(4) Eggen (1969) 


(5) Leckronc (1973) and references therein 

(6) Cowley et al. (1969), Jaschek et al. (1980) 

(7) Johnson (1966), Jaschek et al. (1980) 

(8) HgMn stars (Leckrone 1973) 
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TABLE 2 

EXTINCTION BUMP PARAMETERS 


Name 

V 1 (Mm* 1 ) 

X 0 (A) 

y (pm' 1 ) 

A (pm' 1 ) 

Comments 

HD 62542 

4.74 

2110 

1.29 

3.11 

1,2 

HD 29647 

4.70 

2128 

1.62 

3.35 

3,4 

Average 

4.599 

2174 

0.992 

5.17 

5, 45 star average 

e^Ori 

4.633 

2158 

0.835 

2.43 

5, A.Q* 1 (max) 

HD 38087 

4.559 

2193 

1.000 

6.68 

5, V 1 (min) 

COph 

4.579 

2184 

1.251 

5.71 

5, y (max) 

HD 93028 

4.625 

2162 

0.768 

2.62 

5, y (min) 

30 Dor region 

4.606 

2171 

0.894 

2.62 

6 

LMC average 

4.608 

2170 

0.994 

4.03 

6 

Seaton curve 

4.595 

2176 

1.051 

5.92 

6 

Extinction curve fitting function: 

E(X-V)/ E(B-V) - A, + A 2 X ' 1 + 

[Ar 1 - 

NOTE: for A.’ 1 < 5.9 pm' 1 , A 4 =A 5 =0 

(A.f- 1 ) 2 /X' 1 ] 2 + V 2 

+ A 4 (AT 1 -5.9) 2 + A 5 

( At 1 - 5.9) 3 


(1) Variances are a(X o _1 )=d0.03, a(y )=±0.04 and a(A)=±0.30 and correspond to the observed range of 
parameters found using both B5 V comparison stars from Table 1. 

(2) Coefficients for the fitting function are A^-2.325, A 2 =1.346, A 3 =2.464, A 4 =0.860, A 5 =-0.006. 

(3) Variances are a(A. o ' 1 )=±0.01, c(y )=±0.09 and a(A)=±0.25 and correspond to the observed range of 
parameters found using both HgMn comparison stars from Table 1. 

(4) Coefficients for the fitting function are Aj=0.358, A 2 =0.725, A 3 =3.975, A 4 =0.783, A 5 =-0.168. 

(5) These results from Fitzpatrick and Massa (1986) indicate the average bump parameter values and the 
most extreme deviations in A^" 1 and y among the objects in their sample. Both HD 62542 and 
HD 29647 have A^* 1 and y lying outside the range found by Fitzpatrick and Massa (1986). 

(6) Data from Fitzpatrick (1986). 
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TABLE 3 

SENSITIVITY OF EXTINCTION BUMP POSITIONS TO PARTICLE 

CHARACTERISTICS 


Substance 

size( pm) a 

shape 

X Q Vm' 1 ) 

VA) 

TKpm' 1 ) 

Notes 

graphite 

r<0.005 

sphere 

4.83 

2070 

1.19 

1 

it 

0.01 

tt 

4.76 

2100 

1.19 

1 

ti 

0.02 

ft 

4.54 

2200 

1.26 

1 

ti 

0.03 

it 

4.35 

2300 

1.71 

1 

ft 

0.04 

tt 

4.21 

2375 

2.62 

1 

tt 

"small" 

spheroidal(a/b=2) 

5.10 

1960 

1.19 

1,2 

ti 

it 

" (a/b=1.67) 

5.00 

2000 

1.19 

1,2 

it 

tt 

" (a/b=1.0) 

4.83 

2070 

1.19 

1,2 

if 

tt 

" (a/b=0.5) 

4.54 

2200 

1.19 

1,2 

amorphous carbon 

"small" 

sphere 

4.17 

2400 


3 

amorphous carbon 

annealed at 2200 C 

"small" 

sphere 

-4.55 

-2200 


3 


a For particle size, "small" implies the particle is in the Rayleigh limit . 


(1) Bulk optical constants for graphite have been taken from Taft and Phillips (1965). When optical 
constants from other investigators are used, the bump positions shift by about 100 A (see MRN). 

(2) Spheroidal particles in the small size limit where a and b refer to the major and minor axis 
respectively. The bump parameter values listed are from Wickramasinghe and Nandy (1974) 

(3) These values for the bump position are from Hecht (1986). 
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Figure 1: 


Figure 2: 


Figure 3: 


Figure 4: 


FIGURE CAPTIONS 


Logarithm of the IUE flux (in ergs cm* 2 s* 1 A* 1 ) data of HD 62542. Also plotted are 
the data for the two B5 V comparison stars listed in Table 1. The gaps in the spectra 
correspond to the removal of reseaus and hot pixels from the data. The absolute flux for 
each object can be retrieved by adding a constant, C, to the scale shown on the Figure. 
The constants for 23 Cyg, p Aur, and HD 62542 are C=0, 1.0, and 0 respectively. 


Logarithm of the IUE flux (in ergs cm* 2 s* 1 A' 1 ) data of HD 29647. Also plotted are 
the data for the two HgMn comparison stars listed in Table 1. The gaps in the spectra 
correspond to the removal of reseaus and hot pixels from the data. The absolute flux for 
each object can be retrieved by adding a constant, C, to the scale shown on the Figure. 
The constants for k Cnc, t) Her, and HD 29647 are C=1.5, 2.25, and 0 respectively. 


The normalized extinction curve determined for the line of sight to HD 62542 using 
the B5 V comparison star p Aur. The curve determined using the comparison star 
23 Cyg (B5 V) is essentially identical as are the curves found using a B3 V and B4 V 
(see § Ilia). The open circles represent the average extinction curve of Seaton 
(1979). Also plotted are the three components that result from the fit procedure of 
Fitzpatrick and Massa (1986; see also § IVa) and the fit to the observed extinction. 


The no rmalize d extinction curve determined for the line of sight to HD 29647 using 
the HgMn comparison star u Her (B9 Hip). The curve determined using the HgMn 
comparison star k Cnc (B7/8 Hip) is essentially identical (see § Etta). The open circles 
represent the average extinction curve of Seaton (1979). Also plotted are the three 
components that result from the fit procedure of Fitzpatrick and Massa (1986; see also 
§ IVa) and the fit to the observed extinction. 
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Figure 5: The computed best-fit Drude profiles normalized to a peak value of unity for the 

extinction curves observed towards a) HD 62542 and b) HD 29647. The filled squares 
represent the observed bump profile (binned in 20 A intervals), normalized to unity, 
found by subtracting the computed linear background and FUV extinction fits from 
the curves shown in Figures 3 and 4. Also plotted is the Drude profile for the average 
interstellar curve (Seaton 1979) as well as the extreme range of X 0 _1 found by 
Fitzpatrick and Massa (1986) for 45 lines of sight including a number exhibiting very 
peculiar extinction. 

Figure 6: The environment through which the line of sight to HD 62542 passes. The directions to 

the nearby hot radiation/wind sources £ Pup (05 Iaf) and y Vel (WC8 + 09 1) are also 
shown. The inset at the top shows (to scale) the observed distribution of dust as well as 
the position of HD 62542 (open square) and the distribution of observed IRAS point 
sources. 

Figure 7: Comparison of the fit to the extinction curves observed towards HD 62542 and 

HD 29647 to the fits for the curve for the 30 Dor region, the LMC average curve, and 
the average Milky Way curve (Seaton 1979). These curves were generated from the fit 
parameters found in Fitzpatrick (1986). 
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